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Schwarzschild black hole obey the estimates 101 < CsVj^^^^ and \dt(f>\ < Cgvl^'^'^^ on a compact region 
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1 Introduction 



A major open problem in general relativity is that of the nonlinear stability of Kerr spacetimes. These 
spacetimes are stationary axisymmetric asymptotically flat black hole solutions to the vacuum Einstein 
equations 

Rfif = 

in 3 + f dimensions. They are parametrized by two parameters (M, a), representing respectively the mass and 
the angular momentum of a black hole. It is conjectured that Kerr spacetimes are stable. In the framework 
of the initial value problem, the stability of Kerr would mean that for any solution to the vacuum Einstein 
equations with initial data close to the initial data of a Kerr spacetime, its maximal Cauchy development 
has an exterior region that approaches a nearby, but possibly different, Kerr spacetime. 

Kerr spacetimes have a one-parameter subfamily of spacetimes known as Schwarzschild spacetimes for 
which a = 0. The Schwarzschild metric in the so-called exterior region can be expressed as 

where dugi denotes the standard metric on the unit sphere. In view of the nonlinear problem, it is conjectured 
that a spacetime that is close to Schwarzschild initially will approach a Kerr spacetime that is also close 
to Schwarzschild, i.e., a <ti M. In other words, we can consider the stability of Schwarzschild spacetimes 
within Kerr spacetimes. (Notice that the Schwarzschild family itself is not asymptotically stable since a Kerr 
spacetime with small a can be considered as a small perturbation of a Schwarzschild spacetime.) 

To tackle the nonlinear stability of Schwarzschild spacetimes within the Kerr family, it is important to 
first understand the linear waves 

on the exterior region of Schwarzschild spacetimes. This can be compared with the nonlinear stability of 
Minkowski spacetime whose proof requires a robust understanding of the quantitative decay of the solutions 
to the linear wave equation [5], [19]. 

The pointwise decay of the solutions to the linear wave equation on Schwarzschild background is proved 
in [ini, [4]. In particular, Dafermos-Rodnianski proved a decay rate of < C (max{l,w})~^ everywhere 
in the exterior region, including along the event horizon |10j . The subject of this paper is to improve this 

— -4-6 

decay rate. In particular, we will prove that for arbitrarily small 6 > 0, |0| < Cs^r (max{l,u}) ^ in the 
region {r < R} for any R > 2RI, including along the event horizon. 

Our proof applies a new vector field commutator S that is analogous to the scaling vector field in 
Minkowski spacetime. We will show that for solutions to Ogcj) — 0, Dg (50) decays sufficiently towards 
spatial infinity and only grows mildly towards event horizon. We then prove energy estimates for 5*0 with 
the help of (a slightly modified version of) the energy estimates of <p in [10]. This will enable us to prove the 
decay of Scj). With this decay, we follow Klainerman-Sideris [15] to improve the decay rate for dt(t>- We also 
introduce a novel method to improve the decay rates for (f) and its spatial derivatives. 

We hope that this improved decay will be relevant for nonlinear problems. We recall for example the 
wave map equation from M'^'^ to given by: 

= ((at(/.)' - |V(/.p) . 

To prove the global existence for small data for this equation, it is insufficient to have \d(j)\ < C {1 + \t\)~^ ■ 
One needs an improved decay \d(j)\ < C (1 -|- |t -I- r\)~^ (1 -I- |t — r\)~ . Moreover, one needs the nonlinearity 
to satisfy the so-called null condition (see [13]). In a future work, we will use the improved decay rate we 
prove in this paper and study the global well-posedness of small data for a nonlinear wave equation satisfying 
a null condition on a fixed Schwarzschild background. 

In Section 1.1 and 1.2 we will introduce the Schwarzschild spacetime and the class of solutions that we 
consider. This will introduce the terminologies necessary to state the main theorem in Section 1.3. We will 
motivate our proof with a comparison with the linear waves on Minkowski spacetime (Section 1.4). We then 
mention some known results on linear waves on Schwarzschild spacetime (Section 1.5). We especially discuss 
the work lOJ whose techniques are important for this paper. We will then provide some heuristics for our 
proof of the main theorem in the final subsection of the introduction (Section 1.6). 
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1.1 Schwarzschild Spacetime 

Schwarzschild spacetime is the spherically symmetric asymptotically flat solution to the vacuum Einstein 
equations. The Schwarzschild metric in the exterior region is 

g ^ - il- —j dt^ + il- —j dr^ + r^das2, 

where da§,2 denotes the standard metric on the unit sphere. It is easy to observe from the metric that the 
vector field dt is Killing and it is orthogonal to the hypersurfaces t = constant. Spacetimes with this property 
are called static. It is also manifestly spherically symmetry and therefore has a basis of Killing vector fields 
rii generating the symmetry. Moreover, Schwarzschild spacetimes are asymptotically flat. This means that 
the metric approaches the flat metric as we go to spatial infinity (r —^ oo). 

Synge 23J and Kruskal [17] showed that the Schwarzschild metric can be extended past r — 2M as a 
solution to the vacuum Einstein equations. Its maximal development is usually described by a Penrose dia- 
gram, which depicts a conformal compactiflcation of the 4 dimensional manifold quotiented out by spherical 
symmetry (Figure 1). In this diagram, the coordinate system {t,r > 2M,uj E with the metric described 
above represents the region I, which we will call from now on the exterior region. In the nonlinear stability 
problem, it is this region that is conjectured to be stable. Extended beyond r = 2M, the Schwarzschild 
spacetime contains a black hole (region II in the diagram). Physically, an observer outside the black hole 
region cannot receive signals emitted inside the black hole. The null hypersurface r = 2M separating the 
exterior region I and the black hole is known as the event horizon T-l~^ . 




Figure 1: Schwarzschild spacetime 



We return to the discussion of the exterior region of the Schwarzschild black hole. For notational conve- 
nience, we let 

2M 
^ = . 

r 

We denote as r* the Regge- Wheeler tortoise coordinate 

r* r + 2M log (r - 2M) ~ 3M - 2M log M. 
In these coordinates, the Schwarzschild metric in the exterior region is given by 

5 = - (1 - ^) + (1 - /i) dr*'^ + r^dcrs2 . 

Notice that in the above equation we have used both r* and r. Here, and below, we think of r* as the 

coordinate and r as a function on Q, with r (q) = \J^^^-i i-e. the physical radius of the 2-sphere under 
which the metric is symmetric. The coordinate r* is -|-oo at spatial and null infinity; — oo at the event 
horizon and at r = 3M. The set {r = 3M} is known as the photon sphere. On this set trapping occurs: 
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there exist mill geodesies that lie in this set. In particular, these geodesies neither eross the event horizon nor 
approach null infinity. This suggests, via geometrical optics considerations, that one has to lose derivatives 
while proving energy estimates. We will return to this point when we discuss the vector field X. 
We notice that as in the coordinates (t.r.uj), dt and il are Killing in the {t,r*,uj) coordinates. 
We also define the retarded and advanced Eddington-Finkelstein coordinates u and v by 

t = V + u,r* = V — u. 

At the event horizon H+, u = +00. At future null infinity 1+, v = +00. Notice that in these coordinates, 
the metric is given by 

—4 (1 — /i) dudv + r^da§2 . 
In particular, this shows that 9„, 9„ are null. 



1.2 Wave Equation and Class of Solutions 

We would like to study the solutions to Dgf/) = in the exterior region of Schwarzschild. Written in local 
coordinates, 

ng(f> = - (1 - m)"^ d^(t> + (1 - m)"^ ^--29^. (r'^dr^cp) + ^(p, 

where ^ denotes the Laplace-Beltrami operator on the standard 2-sphere with radius r. 

We notice that Dg commutes with Killing vector fields. In particular, OgCp = implies Dgdtcj) = and 

Dgflcj) = 0. 

The decay result that we prove apply to solutions to the wave equation that is in some energy class initially. 
We define the energy classes using currents of vector fields. We will briefly introduce the relevant concepts 
here in order to present the energy classes. A more detailed description of the vector fields will be presented 
in the next section. 
Define the energy-momentum tensor 

T1J.U = dfj,(l)d„(j) - ^g^ud^^dacj). 
Given a vector field V^, we define the associated current 

Jj (</)) = VT^. (</)) 

and the modified current 

To define the energy classes we need two vector fields: 

i fj, 

Z = v?du + v^dv, 

where yi, 2/2 > are supported near the event horizon with 1/1 = 1, 1/2 = at the event horizon. The precise 
form of yi , j/2 will be defined later. We also define a modifying function for the associated current of Z: 

^ "itr* (1 - 11) 
r 

We note here that 

J {<!>) nldVok, ^ p £ (^^^ + {d^cpf + ly^p j r^dAdr\ 
~ f (^u'{ducpf+v\d,<l>f + {l-n) (^{u' + v') + ^X^^ '^j r'dAdr*, 
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where dVoltg is the volume form of the shoe {t = to} (See Section 2). 

Let S = tdt + r*dr-. 

Define 



Eo{(j> 



.k=0 k=0 J 

1 4-m 

m=0 fc=0 

2 2-m 

m=0 fc=0 

1 1 

^^^;o(a^^^V)+i^l (</>), 

m=0 fc=0 

2 4-m 2 2-m 



m=0 fe=0 



m=0 /s=0 



We notice that the boundedness of these quantities should be thought of as requirements of regularity and 
decay. In the above, Eq, E\, E2, Es and E4 requires 4, 8, 10, 8 and 10 derivatives respectively. In terms 
of spatial decay, all the energy classes require decay of <j) at spatial infinity. However, we note that (j) is not 
required to decay toward the bifurcate sphere {H~^ fl W~ in Figure 1). In the following, we will work as 
if is smooth and supported away from spatial infinity. This assumption can be removed by a standard 
approximation argument. 

1.3 Statement of the Main Theorem 

We prove both pointwise decay and energy decay for solutions of Og(f) = 0. Prom this point onwards, we 
assume > 1, > 1. 

Main Theorem 1. Suppose (j) is a solution to the wave equation on the exterior region of Schwarzschild, 
i.e., Dg(p = 0. Then for any 5 > and any 2M < R< 00 

1. Pointwise Decay of <p 

{v,,u) I < Cs^HV^^^^Ej (0) for r<R. 

2. Pointwise Decay of Derivatives of (j) 

\VcP{v.,u)\<Cs,Rv;^^^ E El{drn''cp) forr<R, 

k+m<l 

where V denotes any derivatives. 

3. Decay of N on- degenerate Energy in the Region r <r2 



r {V<P {t,)f dVolt, + [ (-^ {ducl>f + (1 - m) \f<t>?] d^rfM{.=..} 



fe+m<l 



for any r'l, r'2, r'l. 



5 



Remark 1. The integral in statement 3 represents the part of non- degenerate energy restricted to the region 
r < r2 (See Section 3.5). It should be compared with the corresponding part of the (degenerate) energy 
generated by the vectorfield T = ^ 

r {W<l>(t,)fdVolt, + I ({d^cj^f + (1 - /i) ly^p) dAdU{,=,^y 

Jrl J{r'<ri} ^ ' 

The non- degenerate energy provides strong control of the solution near the event horizon 7i+. Such control 
also easily allows one to extend the obtained decay results past the event horizon. 

For the time derivatives, we have better decay estimates: 

Main Theorem 2. Suppose (p is a solution to the wave equation on the exterior region of Schwarzschild, 
i.e., Og(t> — 0. Then for any d > 



1. Pointwise Decay of dt4> 



\dt^ {v*) I < Csv:'+'eI (0) for r* < I 



2. Decay of Non-degenerate Energy of dt4> in the Region r* < ^ 

{Wdt4>{U)fdVolu + I {-^ {duMf + (1 - /i) \fdt4>\A dAdu{,^,^} 

J{r*<r|} V ^ M / 

<C5^,.min{i,,t-4-4+5^3(0). 
for any rl, r^. 

We would like to point out that the pointwise decay rates in both theorems apply to region of finite r 
and along the event horizon. 

1.4 The Case of Minkowski Spacetime 

At this point, we would like to discuss some decay results for the linear wave equation on Minkowski 
spacetimes. We would like to especially highlight techniques that are relevant to our result. In Minkowski 
space M'^'^, the solutions to the wave equation with initial conditions (p {t = 0, x) = and dt(f> (t — Q,x) — 4>i 
can be written as 

(j) {t, 2^) = ^ (dt J^^ tcpo {x + ty) dA (y) + tc^^ {x + ty) dA (y)^ . (1) 

This formula implies immediately that \<j){t,x) \ < j^, where t — max{t, 1}. This decay is optimal in the 
variable t. However, improved decay can be seen in the null coordinates v = ^ + f^nd u = ^ (t — r), 

3 

where = x'f . In particular, ([1]) implies the strong Huygens' Principle, asserting that (p with compactly 

4=0 

supported initial data is compactly supported in the variable u. Therefore, denoting v+ = max{w,l}, 
= maxjw, 1}, we have in particular 

101 <-^, V7V>0. 
If we just focus on the region {r < |}, where t ^ v ^ u, the decay can be written as 

\<P\<^, ViV>0. 
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However, the use of the representation formula ([T|) is not available on perturbations of the Minkowski 
spacetime. In [5], |19| . a more robust understanding of the decay of the linear waves was necessary. This 
was achieved by the vector field method. Let (f) he a solution to the linear wave equation on Minkowski 
spacetime, □m0 = 0. Define the energy momentum tensor 



Notice that the wave equation implies that the energy momentum tensor is divergence free, i.e. 

V^T^, = 0. 

Given a vector field V^'', we define the associated currents 

tV 



j; (0) = VT^, (0) , 



and the modified currents 



1 ,z„„ ,„ , 1. 



K^'^ (0) = (0) + -w^d'^f^d,^ - -Dgw 



V a2 

9^' 



where is some scalar function associated to the vector field V. Since the energy momentum tensor is 
divergence free, it is easy to check that 

Jj^-' (0) = i^'^'"'" ■ 

Notice that {(f) = whenever V is Killing. In this case (0) is divergence free. Therefore, for any 
solution (j) and Killing vector field V, there is a conservation law: 



Jq ((/))dxt, = / Jo i(j})dxto. 

1 Jt = to 

This is a manifestation of Noether's Theorem, which states that a differentiable one-parameter family of 
symmetries gives rise to a conservation law. We call the vector field V in this application a multiplier 
because we "multiply" it to the energy momentum tensor. An example of this is to take the Killing vector 
field dt and derive the energy conservation law 

{dt^f + E (^-.'^)) dxt = J (^{dtcl^f + Y (d.Af^ dxt,. (2) 

Besides being multipliers, vector fields can also be used as commutators. This means that we commute the 
vector fields with □„!. For example, since d £ {dt,dxi} is Killing, [□,„,(?] = and therefore {d(j)) = 0. 
Then the energy conservation law ^ can be applied to d(f> and we can control the norm of the derivatives 

of (j) of orders 1 and 2. Then using a Sobolev-type inequality ||0||l~(r3) < (which 

holds for compactly supported functions), uniform boundedness of the solutions to the wave equation can 
be proved. The Killing vector fields Qi generating the spherical symmetry can also be used as commutators. 
This is especially useful because compared to the angular derivatives, Qi has an extra factor of r, i.e., f2 ~ r^. 
This allows one to prove in [T^j that for (p decaying sufficiently fast at spatial infinity: 

2 2 
m=l k=0 
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which imphcs a decay in the region {r > |}: 



, ,, C 
< - < — 



after applying ([2]) to fi''^. 

To achieve decay of (/) in {r < |}, one can use the conformally Kilhng vector field Z = {t^ + r^) dt + 2ir9r 
introduced by Morawetz [22]. In this case, {(j)) ^ 0. Nevertheless, by defining = 2i, iC'^''^ ((^) = 
and therefore / J(f (</)) dxt is a conserved quantity. Moreover, some algebraic manipulation would show 

Jo^'"'" (0) dxt > cy (^i;^ {d,<pf + u2 + + + |y0|2^ ^ da;, 

where y denotes the angular derivatives. The conserved nonnegative quantity / J(f {(j)) dxt is known as 
the conformal energy. For the region {r < |}, notice that the boundedness of the conformal energy implies 
a local energy decay 

/ {t^ + {d,cj,f + {ducj^) + \y^Adxt<^. 

J{r<i} V'' J t+ 

After considering the equations (9*^0) — 0, Sobolev embedding would imply the pointwise decay |^| < 

for r < |. Notice that in this region t+ is comparable to v+. Therefore, we have in the whole of Minkowski 
spacetime 

w < £^ 

Klainerman-Sideris TS] showed that more decay can be achieved in the interior region {r < ^} for 
the derivatives of (p. They used the scaling vector field S = tdt + rdr as a commutator. Notice that S is 
conformally Killing and [Dm, S] — 2Dm- In particular, if one has Om(f> ~ 0, then {S(p) = SDm,<j>+2Dm(f> = 
0. Therefore, any decay results that hold for ip also hold for 5*0. Klainerman-Sideris [15] showed that 



^ Wu+dMWmR^) < c (I|5^'/'IIl=(r3) + ||9<^||l=(e3) + ||9'0||l2(k3) + ||9170|U.(r3)) 



E 

By cutting off appropriately and using the local energy decay estimates 



\\9dt^\\L^{{r<±}) < fT in{r-<2} 

since — < -p- in this region. Again, using the Sobolev-type inequality above, one shows that \dt(j)\ < 

in {r < |}. The other derivatives can also be estimated first by elliptic estimates and then the Sobolev 
inequality, since | ju+Sj^^l |i2(-K2-) = | |u+A(/)| 1^2(^2-) by the linear wave equation. Therefore 

in{r<^}. 

We remark that in [T5|, the improved decay in {r < |} can also be proved for the function (p itself by 
inverting the Laplacian. As we proceed to prove the analogous decay on Schwarzschild spacetimes, we will 
avoid doing so. This is because on Schwarzschild spacetimes, it is impossible to invert the Laplacian for 
functions that do not vanish on the bifurcate sphere (7i+ n in Figure 1). 

1.5 Some Known Results on the Wave Equation on Schwarzschild Spacetimes 

We now turn to the corresponding problem for linear waves on Schwarzschild spacetimes. The problem of 
the uniform boundedness of solutions to Ogcj) — on the exterior of Schwarzschild occupied the physics 
community for some time. The first mathematically rigorous result was obtained by Wald 24J for solutions 
vanishing on the bifurcate sphere (7i+ n in Figure 1). Kay- Wald [12^ later removed this restriction and 
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proved the uniform boundedness of a more general class of solutions. They used the energy conservation law 
given by using dt as a multiplier as well as the Killing fields {9t, 0^} as commutators. The decay rates 



|0|<Ci;+\ Vr>2M, 



(3) 



where vj^ = maxjw, 1}, = max{u, 1} and depends only on an appropriate norm of the initial data, 
for sufficiently regular solutions to Og(f> — 0, were proved by Dafermos-Rodnianski in [10]. We note that 
the decay rate ([3]) holds in the entire exterior region of Schwarzschild spacetimes, including along the event 
horizon. In addition to the vector fields in [21],[T2], their approach employed several other (non-Killing!) 
vector fields. One is an analog of the Morawetz vector field Z in Minkowski spacetime. It has an associated 
nonnegative quantity which we will call the conformal energy. It has weights similar to that of the conformal 
energy on Minkowski spacetime so that its boundedness would imply a local energy decay. Another is a vector 
field of the form X = f (?'*) 9, *. The construction of this vector field was motivated by [18 . Unlike other 
multipliers, X is constructed so that K^'^ (0) (instead of J-^'^ {(/>)) can be controlled. This is used to 
estimate some energy quantity integrated over spacetime, in particular error terms from the "conservation 
law" of the conformal energy. The estimates of X are iterated together with that of Z to achieve the 
boundedness of the conformal energy. This then implies the decay of (f> away from the event horizon. The 
estimate associated to X can be thought of as an integrated in time local energy decay. It was extensively 
studied in [T],[2],I1],[7],[TD],[2T]. 

In addition, [TU] introduced a new - red shift vector field - which takes advantage of the geometry of the 
event horizon and is used crucially in proving the decay rate close to and along the event horizon. This vector 
field is one of the few stable features of the Schwarzschild spacetime. In particular, it can be used to give a 
more robust proof of boundedness of the solutions to the linear wave equation on Schwarzschild spacetimes. It 
also plays key roles in the boundedness results for the linear wave equation on small axisymmetric stationary 
perturbations of Schwarzschild spacetimes and in the decay result for the linear wave equation on slowly 
rotating Kerr spacetimes [8], [9]. As we will see later, it will make a crucial appearance in this article to 
achieve the improved decay rate along the event horizon. 

The study of pointwise decay was carried out independently by Blue-Sterbenz [J] . They showed a similar 
quantitative decay result for initial data vanishing on the bifurcate sphere, with a decay rate that is weaker 
than [10_ along the event horizon. In the proof they used analogues of the vector fields Z and X but not the 
vector field Y. Strichartz estimates for solutions of the wave equation on Schwarzschild background were 
shown in [21] . We refer the readers to Sections 3 and 4 in [9] for further references on this problem. 

Considerable attention has also been given to the problem of decay of solutions of the wave equations 
on the Schwarzschild spacetime restricted to a fixed spherical harmonic (jji arising in the decomposition 
(j){t,r,uj) — y <f>i{r,t)Yi{uj), lo G S^. Such results for a fixed spherical harmonic have been obtained in [6], 



[TT] . [16], [20]. We refer the readers to Section 4.6 in [9] for a more detailed discussion. 
1.6 Outline of the Proof 

Our proof uses ideas from Dafermos-Rodnianski [10 and Klainerman-Sideris [15]. In addition to the argu- 
ments used in [TO], we introduce a vector field S = tdt + r*dr* which is analogous to the scaling vector field in 
Minkowski spacetime. Since Schwarzschild spacetimes are asymptotically flat, S is still an "asymptotic con- 
formal symmetry" generating an "asymptotic almost conservation law" . However, the error terms away from 
spacelike infinity are in general large. To see this more concretely, we recall that on Minkowski spacetimes. 
Omit) = implies {S(j)) = 0. This does not hold in Schwarzschild spacetimes. Nevertheless, for Ogcf) = 0, 
we still have a (schematic) equation Dg (Scj)) = h{r) (V(/) -I- V^0) with ft, ^ as r ^ oo. The strategy is 
then to go through the argument in Dafermos-Rodnianski |10| and control the error terms that arise from 
□g {S4)) ^ 0. To do so, we use a slightly modified version of the energy estimates that are available from the 
proof in [IQ]. 




As in later parts of the paper, we define ip = Scj). We would like to prove energy estimates for ip similar 
to those for (j> that are established in [10], except for a loss of an arbitrarily small power of t. A key estimate 
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that will be used to prove the main theorem is: 



t2 



(^2 + (v^)2) ^ {r*)dVol < Cst^'+\ (4) 



where x is some weight and ti < t2 < (1-1) ti- A similar estimate is available with -0 replaced by <f> from |10j 
using the X vector field. In order to prove this, we argue in a similar fashion. We want to show, using the 
vector field X, that for ti, t2 as above 



+ iyipf) X {r*)dVol < Cst^^^^ {coni. energy(V')}, 

where the conformal energy is the current of the vector field Z on the boundary {t — ti}. We then hope to 
show 

{conf. energy (V') at ^2} < C Nconf. energy(V') a.t h} + f / (V'^ + (VV-)^) X (r*) dVol \ . 



We then iterate two inequalities to obtain ^ as in [TU]. 

The main difficulty in actually carrying out the above procedure is that each step is only true modulo 
some error terms that need not be small. These are error terms arising from the fact that ip does not satisfy 
the homogeneous wave equation, but only satisfies an inhomogeneous wave equation, which schematically 
can be thought of as Ogip = h{r*) (V0 + V^(/>). If one applies the vector field method to this equation, one 
would generate an error term of the form 

/ / V^'df,'iph{r*){V(t> + V^4>)dVol, (5) 

1 J -co JS2 

for the vector fields V S {dt,X = f{r)dr-',Z = (^t'^ + (r*)^^ dt + 2tr*d*}. (In practice there is still another 

error term if one uses the modified current, but since it can be controlled similarly, we omit the technicalities 
here.) Applying Cauchy-Schwarz, we can control ([5]) by 

1 1 

/ / [y^fdVol] / / h{r*){{yc^f + {y^<l>f)dVol] (6) 

\Jti J-coJs^ J \Jti J-ocJs^ / 

We control the first factor by some energy quantities of ijj which we are in the process of proving. They 
are set up so that we can estimate them with a bootstrap argument. In order that the bootstrap can close, 
we would need to show that the second factor decays or does not grow as ii,t2 — > 00 (for example with 
t2 = (l.l)ti). The precise rate of decay that is necessary depends on the vector field V under consideration 
and is ultimately dictated by what the bootstrap argument requires. To achieve this, we recall the energy 
estimates derived from the X vector field in lOj. In particular, we have 

{Vcpfx{r*)dVol<C, (7) 



{Wcl>)\{r*)dVol<Ct^\ (8) 

§2 



where x is a weight that decays at spatial infinity. ([5]) gives good control for the second factor in ^ 
for the region {— | < r* < |} as long as h and x behaves appropriately. We will slightly improve the 
weight X from [TUj so that we have, loosely speaking, h{r*) < C (1 + \r* |)~^ x (^*)- This would give control 
for the second factor in ^ for the region { — | < r* < |}. For the regions {r* < — |} and {r* > |}, 
h{r*) < Cil + \r*\)-'xir*) < C(l + t)-' X(r*). Then we can control the second factor in © in this 
region with ([7]) and the extra factor of (1 + t)^^. The reader should keep in mind that these are only 
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heuristics and are not true if directly applied. The actual estimates for these error terms are slightly more 
involved considering firstly that might grow t; and secondly that we do not have energy estimates that 
control every derivatives of '0; and thirdly that some error terms would tend to infinity as r approaches the 
event horizon. The relevant estimates will be proved in Section 5. 

In [ini, the estimates for ip are used to prove the decay for dt(t) in Minkowski spacetime. We show that 
it is possible to argue similarly to prove the decay for dt4> in Schwarzschild spacetimes (Section 7.2). Recall 
that in [15], one then proceeds with elliptic estimates to prove the decay for other derivatives. However, 
on Schwarzschild spacetimes, if we are to prove an elliptic estimate, we are bound to have some lower 
order terms involving only one derivative of 0. These terms cannot be controlled by the estimates of and 
therefore we are unable to use a similar method to prove the decay of the the spatial derivatives of </>. 

Therefore, we introduce in this paper a new method, based on a novel application of 5*, to prove the decay 
for the function cp as well as its derivatives in spatial directions (Section 7.1). We notice that by 



In order to show that the same holds for any t, we note that S is strictly timelike on a compact set of r*. 
Therefore we can integrate in the direction of S from the slice t to a generic slice t. This integration would 
not give an extra factor of t precisely because we already have the estimates for ip = Scj). After controlling 
the spacetime terms by Q and ([5]), we show that for any t, 



We use Sobolev Embedding to get the pointwise decay estimate for <j) and its derivatives (for rl < r* < rj) 
after commuting with an appropriate number of Killing vector fields. We note in particular that in this 
proof, it is unnecessary to invert the Laplacian on Schwarzschild spacetime to prove the decay of (f>. 

The argument above gives the decay of and its derivatives in a compact region of r* , i.e. a compact 
region of space that is also away from the event horizon. (Recall r* is defined so that r* = — oo at the 
event horizon.) In order to prove that also decays along the event horizon, we use the red-shift vector field 
introduced in [10 . This vector field was used in [TU| to show that in some (explicitly identified) neighborhood 
of the event horizon, some energy quantity on an initial slice can control some similar energy quantity in a 
spacetime slab provided that the error terms that are supported in a compact region of r* can be controlled. 
It is then used to propagate the decay of (j) from a compact region of r* to the event horizon. In this article, 
we show along these lines that any decay estimate proved on a compact region of r* can be propagated to 
the event horizon, giving rise to a decay estimate of the same rate. This will be carried out in Section 6 and 
will give the full improved decay result. 

2 Notations 

Before proceeding, we would like to first define the notations used for the coordinates and volume form. 
For the r, r* coordinates, we always use * to denote the Regge- Wheeler tortoise coordinate of the same point. 
For the t coordinates: 

to denotes the time slice on which the initial data is posed. 

denotes the time slice on which we would like to control the solution. 
ti denotes dyadic time slices (which will be defined in Section 4). 
t denotes a generic time slice. 




for ti < t2 < (1.1) ii. Therefore, there exists a time t £ [^1,^2] such that 
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We assume tQ,t^,,ti,t > 0. 
For volume forms: 

dVol denotes the spacetime volume form, dVol = (1 — /i) dAdr*dt. 
dVolt denotes the volume form on a time slice, dVolt = r^^/l — fidAdr* . 
dVoly denotes a volume form on a w slice, dVol^ = \J\ — /ic?Ac?u0 
dA denotes the volume form on the standard sphere of radius 1. 

Whenever we write J without integration limits, it denotes the integration over "whole space" that is ap- 
propriate for the volume form. 

3 Vector Fields 
3.1 Conservation Laws 

We consider the conservation laws for satisfying = 0. Define the energy-momentum tensor 

We note that T^^ is symmetric and the wave equation implies that 

V^T^, = 0. 

Given a vector field V^, we define the associated currents 

(0) = VT^, (0) , 

where tt^^, is the deformation tensor defined by 

In particular, {(j)) = tt^^ = if F is Kilhng. Since the energy-momentum tensor is divergence-free, 
We also define the modified current 

('^) = < (0) + 1 {y^d^<P' - d^w^^) . 

Define K"^-"" {(j)) = (0) + \wd'' ipd^cf) - lUgW<j?. 
Then 

Jj'"' (0) = X^'"' (</.) . 

We integrate by parts with this in a region B bounded to the future by Ei and to the past by Eq. The 
region B should have no other boundary. Denoting the future-directed normal to Eq and Ei by n^o "-Si 
respectively, we have 

Proposition 1. 

/ Jl (0) n^^^dVol^, + f K"" {^) dVol - / (0) n'^^dVol^,. 

JB JT.0 
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Figure 2: Regions of integration 



In this paper, there are two choices of that we will use. The first is to choose to be f = constant 
slices. The second choice is for estimates near the event horizon. In this case, B = {vq < v < vi,t > to}, 
Eo = {u = vo,t > io}U{wo < V < V\,t = to} and Si = {t; = vi,t > to}^{vo < v < v\,u = 00} (See Figure). 

One can similarly define the above quantities for the inhomogeneous wave equation Dgip = F. In this 
case, the energy-momentum is no longer divergence free. Instead, we have 

In this case, 

WJ^ (tP) = (V') + FVd^i;. 

For the modified current, 

V J^''" (V) = K^''^ {ip) - '^Fwip + FVd^i). 

Proposition 2. 

j J^-'" (V') n^/^ofei + j K^'"" W dVol + j [-^Fw^j + FVd^ii^ dVol = j jj'"' (V-) n'^^dVol^,. 

In the case of wave equation on Schwarzschild bac;kground. we can compute the energy momentiim tensor 
explicitly in local coordinates {t,r* ,x^) or equivalently {u,v,x^,x^), where x^,x^ is an orthonormal 
basis on S^. 

Tun (0) = {du<l>f , 

Tuv {4>) = (1 - |y<A|^ 

Taa {<P) + Tbb {<!>) = IM'' - d"cpdc,<p. 

As a result, 

{duV, + - ^ - 14) (|y</.p - 5"<^5„</)) . 

^ Most of the time it is clear from context whether we are integrating over a. t ov v sUce. We will specify in the case of 
possible ambiguity, for example dV ol^t=2v-r*} is the volume form on a t slice, where t has the specifiecl value. 
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3.2 Vector Field Multiplier T 

Define T ^ df Recall that T is Killing. Therefore, 

(0) = 0. 

In the following, we will consider this current on a constant t-slice. 
One computes that in local coordinates 

Jj = ^-^i== [{dtC^f + {dr.cpf + (1 - /i) 1^01^) , 

where is the normal to a t-slice. 



3.3 Vector Field Multiplier X 

Define X — f (r*) dr'. In the following we will use different hmctions /. One computes that 



We consider the modified current using w = 2/' [r*) + ^^^V^/ (r*). Then 



1 + -]f{r*) 

l-ti 2 \ r J 

- i (t^/'" + J/" (ri + (^*) - ^ (3 - V) / in 

where is the normal to a i-slice. 

The vector field X is constructed to control a spacetime integral by the boundary terms, i.e., one hopes 
to control the integral of K^'"^ {(p) by the integral of J^'"' {(j))n^. In order for this to be useful, we 
need {(f) to be everywhere positive. Such vector fields are constructed in [lOj using spherical harmonic 
decomposition. In particular, it was shown in [10 that there exists a family of vector fields Xi — fi (r*) dr* 
for / > such that for any function (f> (not necessarily satisfying the wave equation), if we write out the 

oo 

spherical harmonic decomposition (f) = (f>i, if^''™^' (0;) > 0. 

1=0 

Moreover, one has 



2 a2 



\dA<C K^'^^ ' {(bi)dA (10) 



{{l + \r*\ fil-f,) (l + |r*|)^ (1-^)7 - J,^ 
for ^ > 1, where C can be picked to be independent of I, and 



2 



dA<C I K^'-"" ' (MdA 



(l + |r*|)i+« r2(l-^) 

Moreover for this choice of Xi , the boundary terms are also controllable as shown in [10] : 

K""-^"' {<l,i)dVol <C I Jl{^)nldVok,. 
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Remark 2. We note that although K'^'"'^ ' (0) is shown to be nonnegative everywhere, it has a weight in 
front of that degenerates at r = 3M. Therefore, we cannot directly estimate the integral of I^^P by 

that of K^^'^ ' (0). Instead, we will consider the equation Og {fl(f>) = and estimate the relevant quantities 
with J K-^''"^ ' dVol. This loss of derivative is related to the trapping phenomenon that we mentioned 
in Section 1.1. 

In section 3 we will construct two more vector fields of this form. One will be a modified Xq to control 
a weighted L^-norm of the zeroth spherical harmonic and the other will be used to control the behavior at 
infinity. 



3.4 Vector Field Multiplier Z 

Define Z — u^du + v'^dv This is the analogue of the conformal vector field in Minkowski spacetime. Like the 
case in Minkowski spacetime, it is used to show decay for the solution to the wave equation. One computes 
that 

We consider the modified current using — H£l_lLJi), Then 
2 4r 2r / 8 ^ ^ 



4r 2r / 4 V r 



^ Ay/l — fi \ r T ) 



where n^' is the normal to a t-slice. 



It is shown in lOj that there exist r\,r\ such that for r* < or r* > r-j, K^'^ > 0. 
Moreover, it is shown that / J^''^" n^ dVolt is everywhere non-negative. 
More specifically, if we define S — udu + vdv and 5 = —udu + vdv, 



2 / , \ 2^ 



^ ''m((5</>)' + (^0)') +(1-m) f +(^S^+l<l!^ ] +2{1- fi){u^ +v^)\f(t>\AdVolt 



3.5 Vector Field Multipliers Y', Y and 



Define Y' = ^['lJ du + 2/2 {r*) dy, where yi,y2 > are supported in r < (1.2) rg, with yi = 1, 2/2 = at the 



event horizon and y[ (r*) ^ j/2 (r*) ~ C (1 + |r* |) ^ for 2M < r < r^. 
Here we want to choose tq small enough so that 

1. Y' is supported on r < 3M (i.e. (1.2) < 3M), 

2. K^' ((/)) > on 2M < r < vq, 

3. CK^' U) > -7=JT U) < , , on 2Af < r < ro. 

— ^/l — fj, M ^-^^ {v— const.} — — ^ 

The vector field Y' is designed to capture the red-shift effect at the event horizon [10]. Using the current 
, we will not only produce estimates on constant t-slices, but also on constant v- slices. We will therefore 
record here all the relevant computations. 
We have 



r 
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From this we see that if rp is chosen to be close enough to 2M, requirements 2 and 3 can be satisfied. 
We modify this vector field so that it has better bounds on constant t- slices. 

Define Y ^ Y' +x{r) T, where x ('') is a cutofli' function with X {'"') = \ n. ^ n o^ • ^ ^^"^ following 



properties: 

1. Y is supported on r < (1.2) tq, 

2. K"^ = K^' on r < tq, 

3. CK^ (0) > ^ on 2M < r < r,. 
On the region 2M < r < ?'o, we have 



r>(1.2)ro 



We argue without computation that for ro < r < (1.2) rp, 

|if^i<c^=L=jj(0)4_,j, 

M {t— const.} — ^ ' {t— const.} 

This is true because JT (6) . , controls every derivative of 6 while the terms in jj (ch) . , 

M ' \t — const.} J -r f_i vv/ }v — const.} 

and {(j)) n'^f^^onst } contain only derivatives du and y. Thus the only difference is the weights, which are 
functions of r and are harmless since r is bounded on this region. 

Define N = T + Y . N is clearly causal, thus (0) ^^f^^g^st } — ^' ^^^Y from the horizon, namely when 
r > 1.2ro, {(f>) = (0) "{t=co«st.}- However, as we approach the horizon, j/^ (0) n^^^^^^^^ j - 

J,y (0) n't. . , and thus J,^ (0) n^, , , gives a much stronger bound. Wc assume for our energy classes 

M {i— consi.j {i— const.} o & 

that the integral of ((/>) J^ij^const } bounded initially and this clearly implies the boundedness for the 
corresponding integrals for and initially. The fiux corresponding to should be thought of as a 
non-degenerate energy, which does not degenerate at the event horizon. This allows us to prove decay results 
along the event horizon. 

Before introducing the vector field commutator S, we end this part on vector field multipliers by explicitly 
noting what each of the positive quantities bounds. Most of these are direct consequences of the expressions 
of the currents, except that for J'^'™ , which requires some manipulation and is proved in [lOj . 

Proposition 3. 1. ^ ((a,.0)' + {dtc^f + (1 - j^^p) < CJj (^) <, 

^- I-oo /s^ 7^ (dn^f + (5"<^)' + (1 - A^) + v') \n?) dVok < C J^'-' (0) n>^dVok, 

3- n 7^ ((1 - ^^'^^) dVok < C J,^-^ (<^) <dV-oZ,, 
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3.6 Vector Field Commutator S 

Define S — tdt + r*dr' — vdy + 

This vector field, together with the usual Killing fields, will be commuted with Dg. We note that the vector 
field tdt + fdr is conformally Killing on Minkowski with [D^, ic^t + rd^] — ZD™. Therefore, the commutator 
[□g, S] is expected to approach 20^ towards spatial infinity, where the spacetime approaches Minkowski. 
We set ip — Scf) and derive an equation for ip. 

Proposition 4. 1. [Dg, 5] = (2 + ^i) + f - 1 - ^) a,. + 2 ((^^ - l) - ^) f 



2. DgV = 91 {r*)dr.<l) + 92 (r*) ^(f>, where \gi {r*)\, 
(logr)^ = maxjlogr, 1}. 
Remark 3. Eqmvalently, we write \gi (r*) |, J^^^ ~ (i+|r'|Kiogr)_^ 



I \r*\ r-2M 



Proof. 



[- (1 - ^,)-' dls] = - 2 (1 - ^JL)-' + r*dr. (1 - dl 



[(1 - fi)-' d^,,,S] =2 (1 - fi)-' dl, - r*dr' (1 ~ fir' dl. 



r I ^ r \ r r I ^ \\r I 2r 



2. is immediate from 1. if we let 51 (r*) = 2 (l1 _ 1 _ 2^] and .92 (r*) = 2 ((^ - l) - ^). □ 



4 Estimates for 

The following has been proved in [TU] and is collected for later use. 

Theorem 5 (Dafermos-Rodnianski). 1. J jj^ {4>) rit^dV alt, — J {(j)) Ut^dVoltg, 

2. /^if^-"'"' i^)dVol<C I Jl{4>)nldVolt,, 

I 

3. JJ^{(f>)nldVolt, <CEo{^), 

4- 1% JH W <dVok, < CEo {<!>) t-\ 

2 

5- It^ jK X^'^""""' {(f)) dVol < CEo {(/>) where ii < ^2 < (1-1) h. 
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The following Hardy type inequality is also proved in [10] and will be used throughout this paper. 
Lemma 6. 



{l + \r*\)-^cp\l^li)-^-dVolt„<C / {dr^(lyf{\~ii)—-dVolt, 



Remark 4. This can he written equivalently in local coordinates as 

^r'^dAdr*<c[ [ {dr'(t>fr^dAdr* 

We construct a vector field Xq to control the spacetime integral of (jp itself. 
Proposition 7. 

dVol<C f J^{^)n'^JVolt„. 



Proof. We first notice that we already have control of a weighted L'^-norm of the non-zeroth spherical 
harmonics. This is because by (jlOp . 



/ dA <C [ if^-'^"' i<l,i)dA 

Js2 (l + \r*\t(l-n) - Js- 



{l + \r*\nl-^i) 
for / > 1. This together with Theorem 5.2 would give 

4 dVol <C [ Jl ((j))n'^dVolto 

(l + lr* 1)^(1-^) - J *° 

for / > 1. So it suffices to consider the zeroth spherical harmonic. 
Define Xq = /qS^., with /o (r*) = /„ (r) = -(jq^^ = "sTT+^F^- 
Suppose we act with Xq on the zeroth spherical harmonic of ho- 
using dH), 

i^---^" (^o) .ll^i^ + iiy^oP ( ^) /o irn 
1 - M 2 V / 

^ ^ ^ -/^" (^*) + -/^' (r-*) + 4/o in " ^ (3 - 4/i) /o (r*) ) <I>1 



, . , , /o" (^*) + -/o' in + ^f[> in - ^ (3 - V) /o (r*) , 

1 — M 4\1 — /i r r"^ r'^ / 

^^f"'""" ('/'o) < =^f^/o (^*) 9t0oa..<^o + (fo in + '^^^—^fo in) idtM <t>^, 

where we have used = 0. 

We would have to show firstly that K-^"'"^ ° (0o) > and controls (/)^, and secondly that J^"'"^ ° (0o) is 
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controllable by Jj ((^o) n^^- We first compute the derivatives of fo- 
f'o{r*)={l-ii)drh{r) 

(1+4/^-2)2 - 

= (i_ )^( V 

''^ M(1 + 4/x-2)' (1 + 4/1-2)2 j (1 + 4^-2)2 







_ sf 384 48 ^ , 3/x(l-/x)^ 16r , /xr 

~ (,-'- ~ A*j 



M3(l + 4Ai-2)4 /x(l+4M-2fy r y /x(l + 4/x-2f (l + 4/x-2)2 

/x'(1-m)(3m-2) 
(1 + 4^-2)' 



A computation shows that 



' :/o' + i/o' + ^/o-^(3-4m)/o 



r-'" r2''^ r3 

(192 + At (128 + n (-784 + n (464 + /i (-28 + m (52 + /i (-3 + 4/i))))))) 

4(4 + ^*2)4 

We need to show that 192+/i (128 + /x (-784 + /z (464 + (-28 + (52 + (-3 + 4/x)))))) > for < < 1. 

Case 1: < A* < 1 

192 + 128Ai - 784/x2 + 464/|3 = 16 (-12 - 20/4 + 29Ai2) (/x - 1) > 0. 
52 - 3/i + 4/^2 reaches its minimum at |. Hence, 52 - Sfi + ^jj? > 
-28 + M (52 - 3m + 4/i2) > -28 + > 0. 

Case 2: < At < 

464 — 28/i+ has negative discriminant, hence > 0. 

Also, for this range of /i, 192 + 128/i - > 0. 

Therefore, K^o,w''o ^^^^ ^ ^ Moreover, < CK^°''^''° ((/)o). 

It now remains only to control the boundary terms. Using Lemma 6 and Cauchy-Schwarz, 

J^"'^^" {<Po)n''dVolt 

J \ (1 + |r*|)^ J 

<C [ — ^((a,0o)' + {dr'<Po f) dVok, 

<C J Jl{(l>a)nldVolt.. 



to- 



□ 



We would like to construct a vector field X = f (r*) dr* so as to improve the weights in r of the spacetime 
integral that can be controlled. More precisely, we have the following: 



Proposition 8. 



t^dVolta, 
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to Jl JS^ fe^O-^ 

/or < 5 < i. 

Remark 5. T/ie Zoss o/ derivative above is unnecessary because we are considering only a suhregion of 
{r* > 0}. One can construct yet another variant of the vector field X to achieve the above estimate without 
any loss of derivatives. However, since this would not improve the regularity in our final result, it is not 
pursued here. 

Proof. Let X = f (r*) 9, * , where / = X ('"*) (^1 ^ {i+r')^ ) ^^'^ X is a cutoff function satisfying 



X 

We recall 



r* < 1 

1 r* > max{100, lOOM} 



1 - /i 2 V ^ / 

- I (j^r (rl + -/" {r*) + 4/' (r*) - ^ (3 - 4/.) f (r*) 
4 \ 1 — /i r r-^ 

Since we already have control of the spacetime integrals on a compact set using Theorem 5 and Proposition 
7, we only have to show that i^^.™"" (<^) > for r* > max{100, lOOM}. For r* > max{100, lOOM}, we have 

(1 + r*) 

~ 5a+^H2+j) 

Clearly, the coefficient of (dr-(i)f and in K^^"""^ {(f) is positive for r* > max{100, lOOAf}. We now 

study the coefficient of 0^ in K^'"""" (0) for r* > max{100, lOOM}: 



1 — /i r r'' 

1 (5 (1 + 5) (2 + 5) 4(5(1 + (5) 2MS 



< 



l-fi (l + r*f+^ r{l + r*f+^ r^il + r*)^+^ 

12M 32M2 2n , 
J + J - -T 3 - 4^4 

3(5 (1 + 5) (2 + (5) 4(5(1 + 5) 2M(5 12M 321/^ 



2(i + r*)3+'5 r(l + r*)^+'' r2(l + r*)^+'' r3(l + r*)'' r5(l + r*)'' 



<^(i±4.r^-iU^^r25-i2' 



r(l + r*)'+^ V 2 ) r3(l + r*)* V 100 
<0. 

Hence K^^^"" (</>) > for r* > max{100, lOOM}. 

Moreover, on this region of r* , (r'^'^ idr-(f>f + r-^'^cj)^ + r^^\f(i)\'^^ < CK^''"'' {(j)). 

Finally, we have / J^'"" {(f) n^'dVolt < CJ^ {4>) nf^dVok using Lemma 6 and Cauchy-Schwarz exactly as in 
Proposition 7. □ 
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Remark 6. The weights in the Proposition are the same as those for Minkowski space. Since Schwarzschild 
is asymptotically flat, they are the expected weights. 

Corollary 9. In local coordinates, Theorem 5, Proposition 7 and 8 imply via Proposition 3 the following 
bounds: 

^- ((^'■•'^)' + + (1 - m) IWP) dVok, < C / Jj (0) n^t^dVok,, 

2- I 7^ i9u<l>f + {d,<jyf + (1 - Ai) I WP) dVok, < CE^ (0), 
3. J^T^irn^^2^y^i^^^cEo{q>)t-^ 

I It^ jl^Xnt-l) + (fw^^°^ - ""^'"^^ ti<h< (1.1) ti, 

5- ftl Ih^ (TTFW^^^°^ ^ ^Sto Eo {^''<t>) t-{\ where t, < t, < (1.1) ti. 

5 Estimates for i/j 

In this section, we would like to imitate [TtJ and prove an analogue of Theorem 5. For technical reasons, 
however, we will need to lose an arbitrarily small power of t. 

Theorem 10. 1. J Jj (tA) nldVok, < C {J {^) n'^JVoko + {(j))), 

2. j J^{^)nldVolt, <CttE,{cb), 

3- I[^u^r'<i^}J^i^XdVok,<CE,{c^)t-^+', 

4- Sll jit K^"'"''' (Vj) dVol < CEi (0) t^^+^, where ti < ta < (1.1) ti, 

The general strategy is as follows. We follow the argument in [lOj but now in the conservation law 
for each of the vector fields, there is an extra error term which is a spacetime integral that looks like 
Ito Vi'd^ijngijdVol (as wen as an extra term —jj^' wipOgipdVol for the modified currents). Very often, 
we need to show that this integral decays (or does not grow) with thus we need to "produce" some decay 
in t. We do this by splitting the domain of integration into three regions and estimating them separately: 

1. For the region {| < r* < oo}, we use the fact that Dgtp contains negative powers of r* , (which is a 
consequence of the asymptotic flatness of Schwarzschild) . In this region, negative powers of r* can be 
estimated by negative powers of t. 

2. For the region {— | < r* < |}, we note that we have decay in the spacetime integral of 4> for each 
dyadic slab by Corollary 9.4, 9.5. We therefore estimate the integral on this region by that of {(f>). 
Here, it is essential that we use the improved X estimates given by Proposition 8. 

3. For the region {— oo < r* < — |}, we make use of the fact that there is an extra factor of (1 — /i)^ in 
the spacetime volume form compared to the volume form on a time-slice (see Section 1.5). From the 

definition of r*, we have (1 — /i) < Ce'^'" , thus the factor of (1 — /z)^ gives exponential decay in r*, 
which translates to exponential decay in t in this region. Therefore, on this region, we first estimate 
on each time slice, and then carry out the integration in t. 

Since we will often perform integration dyadically we first set up the notation. We define a dyadic partition 
of [fo)^*] by to < ti < ... < tn — i*, where < (1.1) i^-i and n is the minimal integer such that this can be 
done. In particular, log (t^, — to) ~ n. 
We begin with the T estimate. 
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Proposition 11. 

2 

/ Jj (V) nldVok, <C I Jl (V) n^tJVok, + CY,Eo (^^'<^) • 

■' k=0 

Proof. The conservation law gives 

J jJ(V'X.dFo/t. = J Jl{il,)nt^dVolt,+ j dt^Ug^dVol. 

We split the error term into three parts and estimate them separately. 
By Corollary 9.1, 

/•OO 

I / dtipngXl;dVol\ 



<C sup (r [ J^{i;)n^dVoltY r [ {n'^) n'^JVoh,] f' rUt 

to<t<t, \J-ooJs'^ / \^^qJ-ooJs'^ ) Jta 

1 1 ~ 

<C sup [f I j]^{^)n^dVohY [Y. r I Jli^'^'fXdyoltA 



Vfc=0 

For the middle region, we observe that By Corollary 9.4 and 9.5, 



dt'>l)Ug',\)dV ol\ 

sup jf^ Jj (V) nfrfFo/*) ' I E E ^ (/ ' ^ ~ '^^'''*) 



<C sup [f I Jl{^)n^,dVolX WY^JxT^ r Y^""'^""^^ 
<C sup f /~ / JjWnfdFo/iV (E^o(0'='^) ) 

to<t<U \J -oo JS^ / VrTn / 
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By Corollary 9.1, 

_ ± 

I / ' dt^ag^dVoi\ 

J — oo 

<C sup / j'/^Wn^dVohY r I Jj (i^VXc^^o/*! r e-'^'dt 

to<t<t, \J-oo JS^ / \^J-CX)JS2 J J to 



<C sup 

to<t<t 



These together show that 



j Jl{i^)nldVolt, < J Jl{i^)nt^dVolt,+C ^ sup^ {^j {'il^) n't dV ol^^ (E^oC^'-^) 



which implies the Proposition with the following Lemma, taking h\ {t) = and /12 {t) — ^ Eq (O'' 0) . □ 
Lemma 12. Suppose f {t) is continuous, hi (t), /12 {t) are increasing and we have 



fe=0 



1 , 1 



f{t*)<C{f{to) + hi{U)+ sup fit)^h2{U 

V to<t<tt 

for all >to. 
Then 

f{U)<C{f{to) + hi {U) + h2 iU)). 
Proof. Suppose sup / (t) is achieved by / (t) for some to <i <U. Then 

to<t<t, 

f{t)<C (/(to) + {t}+f{i)h2 (t)^) . 

hi (t), /i2 (t) increasing implies, 

f{i)<C (/ (to) + hi {t,) + f {t)-'h2 (t*)^) . 
Using Cauchy-Schwarz and subtracting \f (t) from both sides, 

f{t)<C{f{to) + hi (t,) + /l2(t*)). 

Clearly, / (t*) < sup f {t) = f {t). Hence we have the lemma. □ 

ta<t<t. 

We then derive an X estimate. Here unlike in the case for 0, in which the X estimate was used to 
improve the already known estimate from Xi, we need to consider both of them at the same time. 

Proposition 13. 

j I ^ (V^) I + ^ ^' (V'i ) dVol 

Jto I 

<c(jjl (V) n^dVok, + j Jl (V') ntJVoko) + Ct^^+' ^0 ("'<^) 



23 



Remark 7. The reader may ask why this Proposition gives decay for the error term, while the statement of 
Proposition 11 does not. In fact, the proof of Proposition 11 is sufficient to show that the error term decays. 
However, we do not pursue this as it is unnecessary for later use. 

Proof. Decompose ^ — \pi into spherical harmonics. 

I 

Since Schwarzschild spacetimes are spherically symmetric, DgV-"; = .9i (^*) dr*<t>i + 52 {f*) f {^(pi). 

Notice that K^'^ [ij)) is not everywhere positive. It is identically zero for r* < 1 and as we have shown 

in the proof of Proposition 8, K^''^ > for r* > max{100,100M}. On the remaining (not necessarily 

positive) region 1 < r* < max{100, lOOM}, we have \K^''^'' (tpi) \ < CK^'''^''' {i^i). (Notice that we have 
avoided the region around r = 3M where this inequality is potentially problematic.) 
In particular, applying Proposition 2 for the vector field X, we have 

J \K^^'"'' {il)i)\dVol + J K^''^''' {tlji)dVol 
< j K^''"'' {7pi)dVol + {C+l) J K^^'-^""' {iPi)dVol 
= j J^'""^ mnldVok, - j {i,i)nldVoh, 

+ {C+1) (^j J^--"" {xPi)nldVok, - j J^^--"' {^i)n^,^dVok)j 

+ \j [f' + ^^^^/) i^inii>idVoi - j (/a,.V;) n,j;idVoi 

+ + 1) Q y [I'l + ^^^^/;) i^i^i^idvoi - j ifidr^i^i) aiJidVo^j 

<C (^J Jj (V;) nldVolt, + J Jj (V;) n'^JVoh, + J | (r" V; + ^r-VO ai,i\dVo?j . 

We split the last term into three integrals and estimate them separately. 
By Theorem 5.2, 

r r I \{r-^i^i + dr^i}i)Ui,i\dVol 

Jto J| JS^ 

<C ' J^J~^^^ (logr)+ {r-i--*\iJi\ + r-i--^\dr*i;i\^ (r'^-J (|9,.(/.| + \fn^\)) dVol 

1 1 

<ct-'^- ( /"' \K^^-' (^0 \dvo^ ' fx: /'' r L 1^^'"'* ("''^') i'^^^'^ ' 



to Ji Js^ I \i.^cJto Ji Js^ 



<Ct-'^- ( j ' r [ \K^'-' (^0 \dVol) IE/' r I l^^''"* (^''^') l'^^'^^ 

\Jto J-ooJs^ J Xj^^Jta J-ooJs^ J 

< / / / 1^^.-* (^,) + Cto E / \ -^l ("''^') <dVoU,. 
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By Theorem 5.5, 
ft, /•$■ f 

\{r-^i^i+dr'i^i)Uiji\dVol 




dVol 



(lif*-'"" (V-i) I {iPi)) dVol\ [y, r j\ I 1^*'"* (^^VO \dVol 




i=0 fe=0 

1 

(lif^^'"" (V'O I + i^^-'""' (V-i)) dFoZ + Ct-^ J2 ^0 i^'^i) ■ 



fc=0 



By Theorem 5.1, Proposition 11 and Lemma 6, 



/ * / V I {r-'^i+dr^-^i)n^i\dVol 

Jto J-OO 

t — - 

<C r ( ' ( \r*\{\lPi\ + \dr^^i\){\dr,C^l\ + \Wl\)dVol 
Jto J-oo JS'^ 



X {r*f (1 - /i) + (1 - A*) I W/P) j 

<C sup f/ Jl{i,i)ntdVolX I Jl{n''(pi)n^tdVok\ e 

to<t<t, / \,._f.J-oo I 



2 



Subtract the terms with K from both sides and get 

^V"^'"" {i^i)dVoi 



Sum over Z > to get the Proposition. 

We locahze the estimates in the above Proposition to obtain decay as in [TO] . 
Proposition 14. Let tQ < h < (1.1) io, \rl\ + jral < ^. Then 



f \K^'^" Vi^^"'""' {^i)dVol 

Jv I 

<C f (V-) n[;dVok„ + to E E ^0 i^^^' 

\ <J 1 n n 
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where V = {h < t < < r* < ^} or V = {to < t < ti,r* - {ti - t) < r* < r* + {h -t)}. 

Proof. Let X = I Q j^ljj > 1 1 ■ Out = to,lct (p ^ X (ofe) ^t't> = X [oMt^) ^t<l> and solve for Dg4> = 
for t > to. 

Following [TU], we have 

^0.715*0 1 I- 

-j==cj>''dVok,< / J^{ct^)nldVolt,. 

-0.715*0 V J- ~ M J 

This is true because of Proposition 3.2, 3.3 and an elementary one-dimensional estimate: 

1/ [x] \Hx < Ca^ ( r \d.f (x) p + 1/ (x) \'dx) , 



1 



for a > lH 

Using this, we can estimate the current of (j): 



Jj (0) n^yVok„ 

0.715*0 /.0. 715*0 -1 



-0.715*0 ' J-0. 715*0 \/l ^ M 

Similarly, 

jJ (r!<^) <dFo?t„ < Cto ^ / {ncf>) nldVolt, 



Ho 



Define 'tp — S4>. By Proposition 13, 



The left hand side equals /p lif'^''"^ (V') I + E {ipi)dVol by finite speed of propagation. 

/ "'t'o "^^"'to '^^'^ tie estimated in a similar way as / ri^^dV^oZtp. More specifically, we claim 



that 



To see this, we first note that 



V' = X ( 777^ ) V' + » X' 



0.65io/ 0.65io V0-65to 



^One can prove this one-dimensional estimate by first considering g = on [— ^, ^] and the trivial bound /"^ M^:^ < 

Ca J^^ (x) Then one sets g (x) = f (x)^ and use Cauchy-Schwarz to get J^^\f (x) \^dx < Ca^ J^^\dxf {x) \^dx.. 
Finally, one cuts off / (x) to be identically zero in [— i, i]. 
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Also note that on the support of ^, ^ < C. Therefore, 

M.7lhta M.ribtQ 1 <.0. 715*0 

< / Jl W nldVok, + Cto' / -^=ij^dVolt, +C jJ (</.) nf/FoZ. 

J-0. 715*0 J-0.715to V ~ J-0.715to 



^0. 715*0 



1 



dVol 



/U. ( 
-y=^cP^dVok, 
-0.715*0 V J- ~ M 

<Cto' (/ W<,dVok, + J (0) < 
<Cto ' (1 W nldVoh, + Eo (</.)) . 

We would now want to control / Jj (^7p^ n^^ dVolt^ ■ Using the conservation law for T and an integration by 
parts in t, 

J Jj(vi) n'i^dVok, 
= j Jl (^) nldVolt, - j dti^Og^dVol 

<jjl (f) n'^jvok, + \ji>ng (dti) dVoi\ + 1 j i>ngi,^fT^dVoit,\ 

+ 1 j i)Ugi;./l^dV0lt,\. 
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We first estimate the spacetime error term in this expression. Using Proposition 3, 7, and 8, 

I J ijOg (dt^') dVol\ 
< J l^pOg (5 (dt4>)) \dVol + J \i;Dg (^dt4>) \dVol 



<C 



I 



(l + |r*|)(logr)+ 



\^\ (\dr.dt^\ + \f(dtnpi \)dVoi 



<c 



J (1- 



{l + \r*\r 



rdVol 



^ (i + \r*[f({dtdr*4^\\f(nd,^A ^ 

/ zr^t —dVol 



74- ± 



<Ct 



k=0 



<Ctl J \K^^-' (^) I + ^i^^' (^,) dVol^ ' 

X ( ^ / {dtn"^'^ I + (dtn'^^i) dvoi 

I iyj Jl {i) <,dVok, + j Jl (V') n^dVok, + to ^0 (^'^) 



k=0 



1 2 



TO=0 fe=0 



where at the third to last step we again used Proposition 13. 

We estimate the boundary terms using Lemma 6 and Corollary 9.1 

I j ^i>Dgi,^^- iidVolt,] 

<C J (l + |r*|)r-3 (logr)+v^r^|Vi| (\dr4\ + dVok, 

<C j Jj (v:) n>i^dVok, + cj (jj (<^) + Jj {n4>)) n^.dVok,. 
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The term for t = ti is done analogously, but with a more careful choice of constant. 

I J i!Dgi!^/l - ndVolt, I 

(1 + |r*|)r-3 (logr)^ VWIV^I {\dr4\ + dVoh, 

Jl (V^) <dVok, +cj (Jl (^) + Jj (f2^)) nldVok, 
=\ J Jl W <dVolt, +C j (Jl (f) + Jj {npi) nldVok,. 
Combining these estimates and subtracting 5 / Jj n^^dVolt^ on both sides, we get 

'' m=0 fe=0 

1 2 

It remains to control ^ ^ -Eg (^il'^d^'tfj . 

m=0 fc=0 
1 2 

m=0 fc=0 

/■0.715to / 1 5 . . ^ / N \ 

/ E E -^.^ (^rf^'^^) < + E E J" {9rn>^^) <„ dv-oio=*c} 

^-0.715*0 \m=Ofe=0 m=ak=0 ) 

/.0.715to / 1 5 14 \ 

/ E E ^ (^™"' >) + E E ^ (^™"'?^) < dvoi^,^,,^ 

J-0.715to \„j^ofe=0 m=Ofc=0 / 

/.0.715to / 1 5 14 \ 

+ Cto' / EE (Sr^^V) + E E^o (Or^'^f r^dAdr* 

J-O.nhto \TO=Ofe=0 TO=Ofe=0 / 

1 2 

^^EE^o(9rfiV)- 



TO=0 fe=0 



After establishing the X estimates, we turn to the Z estimates for tf}. 
Proposition 15. 

j J^^'-" {i>)nldVok, 
<C [ J^^'^\^P)nldVok,+cY, /'* I"' tJ^K''"'"'" {n%)dVol 

J u—n J to J r,* r 



+ c i^fy j\ e^^' {\K^^-'' (V) I + x^i^^--^' (Vz) j dT/o/j ' (e E ^0 

2 



fe=0 
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Proof. By Proposition 2 applied to the vector field Z, 

= I J^'""' (V)</yoZ{t=t,} + j K^'""" {il;)dVol- I ^llS^ZJ^O,i,dVol 
+ j (u^dutp + v^dytp) DgtpdVol. 

As remarked before, there exists rj,r2 with < r| such that K^''^^ (tp) is non-positive for r* < or 
r* > Tj. Therefore, 

J K^'""^ {ip) dVol < ' K^'""^ {ip) dVol 

<C r /" % (V^ + lyVP) dVol 



'to Jr'i 

1 rt. 



For the first error term, we again estimate by looking at three separate regions. 
By Proposition 3.3 and Corollary 9.1, 

-^t„<t<t (/'^/f'"' (V'X'^^o^*)' (e/ jj (fiV)<o'^^««{*=to} j 
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By Proposition 3.3, Corollary 9.4 and Corollary 9.5 

ft. 

I to 



<C sup ( / Jf''-' (i/>) nfdVoi, 



<C sup f /j^^'-'Wnfrfyo/tV fxiE*r'i^o(f2V)l 

to<t<t,\j J Kr^tr'a / 



,i=0 fe=0 
2 



By Proposition 3.3 and Corollary 9.1 and using the fact that (1 — /u) < Ce'^ 



to oo 



,2 



<C sup (^J J^^"""" (V'XdFo/t^ ' ^ * 6-^=* (^E/ -^J ("'''^) 

<c ^ sup^ (^y J^'-' (VXdyoZt) ' {^En (f^V)^ • 

The estimation of the second error term is slightly more involved because there is a factor of in the 
integrand. In particular, even near spacelike infinity, one needs to use estimates for the spacetime integral 

for (j). We intend to estimate this term separately in three regions as above. However, for technical reasons, 
we will divide the regions slightly differently. Divide as usual the interval into < ti < ■ ■ ■ < tn = t*. 

We then set the three regions to be \ \{ti < t < ti+\,r* > -^}, 1) {ij < t < tj+i,— < r* < -^}, 

1=0 i=0 

U{ti<f<t,+i,r* <-^}. 
i=o ^ 
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In the region [J {ij < t < ti+i,r* > -^}, we estimate one power of t by that in J^'"" (0) and the other is 

4 = 

canceled with the decay in r. To achieve this we use Proposition 3.3, Theorem 5.1, 5.2 and Proposition 8, 

n-l .ti+i 



J2 / /. {u^du^ + v'^d^^)agi}dVoi\ 



<Ci|^ sup (^j J^'^"" {ilj)n^dVolty Jj ((/.X/T/oZ: 



For the region {ti < t < ti+i, — -^^ < r* < -;^}, we first rewrite into (t, r*)-coordinates and then perform 

1=0 

an integration by parts in t. It is to avoid extra boundary terms during this integration by parts that we 
have divided our regions differently from before. The reason that we perform this integration by parts is that 
instead of a spacetime integral term with i9t^, we would prefer a term with ^, which can then be controlled 
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by the integral of \K^''^'' \+Y^K^' 

I 

/-ti+i 



<c\^ [ r I \tr*dr>i}Ug^l)\r'^ {I- n)dAdr*dt 

n— 1 t- — 

+ ['^' []. I (f + {r*f)dti'ngi,r'^{l- ii)dAdr*dt\ 

n—l t- — 

r I \tr*dr*i)Ugi}\r'^{l- n)dAdr*dt 
+ r ' r [ t^li^ng {dti)) (1 - /u) dAdr*dt + C V / V / ^l^^s^k^ (1 - /«) d^rfr-*dt 

+ / t^\i'Ogi,\^/T^dVolt, + / ^ t>ngVlv/i^t^^o/t,+, 

j=0 i=0 "'""2- 

<C ' I [ \tr*dr->^Dgij\r'^ {1- n)dAdr*dt + C [ ' [ t^lipOg {dt^)\r^ {1 - dAdr*dt 




n ii 

t\^l)Dg^l)\r^{l- n)dAdr*dt + cY^ /" ' t^i\^l>agi;\^/T^dVolt,. 

i=o -^-^ 

We now group this into three parts: firstly, the spacetime term that grows like t^; secondly, the spacetime 

terms that grow like t; and finally, the boundary terms. 

By Proposition 7, 8, Proposition 3.4, 3.5, Theorem 5.1 and 5.2, 

f r [ t'^\iljDg{dtip)\r^ {1- n)dAdr*dt 

J to J-^Js^ 

p t^logr)^{l + \r*\)\i;\{\dr'<l>t\ + \fm) 



'to 



<c I r r w^^^^dvoi 



Ito J-i {l + \r*\f 



<c{r /' w(|jr*."*wi + j;jr^'-«'"(«)<jra)7|;(-') {j2Y^E„(dTaS 



; / / \i=0 / \m=Ofe=0 



<c[f' 1"/+^' (iX*--^ I H-E^'"' dVoX [j^jZEoidT^^ 

\Jto J-i \ I J J \m=Ofc=0 
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By Proposition 7, 8, Proposition 3.4, 3.5, Corollary 9.4 and 9.5, 

/ / \tr*dr*tpDgtp\r'^ {1- i2)dAdr*dt+ [ [ t\ipDgip\r^ {1 - fx) dAdr* dt 



'to 



I to 



ti+i 



U-^ + \r*\f({dr'^f + \fn<p\A y\ 



[S'- U'v: ^ 



'to 



Kk=0 / \i=0 



By Proposition 3.1, 3.3 and Theorem 5.4, 

n „ li 

i=o •'-^ 

<C^ sup (^J J^^-' (V) nfdFo^t) ' (j^ Eo (f^V)) l) 



fe=0 / \i=0 

i / 1 



<Ctl sup f / J^-""" (^) n^dyo/t) ' ( E ^0 
to<t<t, \J J y^^Q 



These together give: 



\ [ ' {u'^duip + v'^dyilj)n-iljdVol\ 
Jto J-i 

<c I f' ^2+2^ j li^^.™" I + ^i^^'.-"' (^;) j dvoi\ ' E (^r^^ 

\ ; / / m=Ofc=0 

+ ctl ^ sup^ J^'-" (V) <dyoZt^ ' f E ^0 ("''^) ] • 



n-1 



We finally look at the third region, \^{ti <t < fi+i,r* < —-^}-, for the second error term. By Proposition 



i=0 
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3.1, 3.3 and Theorem 5.1. 

t —- 

\f f ' {u^du^ + v^dyij)aijdVol\ 

J to J — oo 



to J — oo 

X ' {r*f (1 - /x)^ {{dr'4>f + I Wl') dVol}j dt 
<C f' e-=* if ^ (^) n^VolA [Y,f ' (^^V) ntdVolA dt 

J to \J-00 J \f,^QJ-00 ) 

<C sup ( [ J^'-" (V) n'^dVolt] ' {Eo (<^))^ . 



Therefore, 



1 j^'-" 

+ c[f' 1"/+^' [\K^'-' m + Y^K^--''' (V'o) dFoA7^^^;o(5^^2^ 

1 

— / 2 \ 2 

sup tl (^J J^'""^ {i))n'^dVolty \^^Eo{fl''^)j . 



+ C 

to<t<t 



The proof concludes with Lemma 12, taking 



2^1 2 



\fc=0 / 

We notice that hi (t) and /i2 (t) are increasing. □ 

We now combine Propositions 11, 13, 14 and 15 to prove Theorem 10.2. This will then imply the other 
parts of Theorem 10. 

Proposition 16. 
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Proof. We first show that / J^'"" {ip) n^,dVolt, grows only hke t\'^^ . Using Propositions 13 and 11, 

/** (V') I + ^ if""""" (V'i) dVol 

Jto I 

I — n 



fe=0 

Similarly, 

r** 

I " /c=0 

Apply Proposition 15 to get 



1 rt> 



2 

1 3 

<C I J^'-" (V) -^l (^'^) ^o'^^o^^o + Ci* E ^0 (^''^) 

■' k=0'^ k=0 

+cti+' ( f j^{^)n^jvokX f E E^o {dr^'A ' +cti+' E E^o (^r^ 

/ P ( ^ C ^ 3— m 

<c / j^^'"'' (v-) <dv^< + E / < ("'^) K'^^^^io + E E ^0 (^r"'<^) 



We now have some control over / J^'™^ (ip) n^^dVolt, and we will use Proposition 14 to estimate the space- 
time integral terms by integrating dyadically. 
By Proposition 14 and 11, 



/ / Y.^^" {^i)dVol 

\ ^ \m=0 fc=0 / / 

<C [tf j 4 Wn^^dVok,+t-'+' / Jj {n'i,) nldVok, + E E [dT^^'A] , 

\ \fe=0 m=0 fe=0 / / 

where here we have not kept track of the constant factor in front of 5, but just note that it can be chosen to 
be arbitrarily small. 
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We can apply the same argument to J^''^^ JJ? ^^ii'-^i'"' ' [flipi) dVol to get 

^ I 

/ / ^K^''"^ ' {fli^OdVol 

\ •' \/c=0 m=0 fe=0 / 

This in turn provides more control on / J^'^^ {ip) n^^dVolt, by Proposition 15: 
j J^^'-" {i,)nldVok, 

fc=0 i=0 •^''1 ( 

2 

+ C4E^o(^'''^) 

fe=0 

1 » /I 4-jn \ 

<^E / -^/f'"' (""v-) n'^jvok, + cti £0 (V) + E E ^0 (^r^^V) 

fe=0 V m=0 fc=0 / 

n-1 I . \ 5 / 1 2 \ 5 

+ ^ E / / , 1^^""' (^) I + E (V'O rf^^o^ E E ^0 (^r^''^) 

i=0 \-^*» -^^1 I / Vm=Ofe=0 / 

«-l / + ! „| . \ 3 / 1 2 

<ci;i (</.) if + c ^ ti+<^- / / ^ li^^'-" (V) I + E (^0 E E ^0 (9ro 

1 4-m 

Here, we recall that we have defined Ei (0) = + ^ ^ Ef) (9^17^) in Section 1.3. 

Clearly, we can replace 5 hy e with a different constant C which depends only on e: 

<CE^{<l>)tl 

n-l I ^t,+i „ I . \ I / 1 2 \ 5 

+ ^ E / / , i^*-'"^ (V') I + E (V'O civ^o^ E E ^0 (^r^'^'/') 

i=0 V"'** "'-I i / \m=0/c=0 / 

Notice that at this point, the only term that exhibits more growth than expected is 

/ ch + i - - \ ^ / 1 2 



n-l / + i ,.L ^ \ 2 / 1 2 



We will close the argument with a bootstrap. 
For notational purposes, we define 



t ■ — 

lit, = / V 1^^'"* (v-) I + E^"""""' (^') 
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prjl is equivalent to 

/*. < C Ue, {4>) + 2 4+' {IhS^ El (<^) ) ■ (12) 



\ i=0 

On the other hand, Proposition 14 gives 



IIu < C (i-'/t. + t-'+'Si (0)) . (13) 



Assume It < At^Ei (0), where A > AC. We want to show that It < ^t^Ei (0) for all t > 
^From the assumption and (jl3p we have 

< C {Atf+'Ei (0) + t-^+'E, (0)) . 

Hence, by picking e = 4 in (fT2l) . 



40oci f A^+l") 



It <CnEi {4>) + lOOC^ [a^ + l) iTi;^ (0) 
iCt^Ei (0) + lOOC^ + l) t^Ei ((/)) 



<( 

<^t'E, (0), 

4 

since A > 4C and i > \ ^ . □ 

Remark 8. We would like to note that the number of derivatives used in the above argument is highly 
wasteful (we used a total of 8 derivatives!). Blue-Soffer constructed a vector field to control trapping with 
only e derivatives J^. Therefore, we can, in principle, repeat the above argument noting the unnecessary loss 
of derivatives. The details, however, have not been pursued. It is known that with this vector field. Theorem 

5 holds with Eq {(f)) only having 1 + e derivatives. Moreover, in Proposition 11-15, instead of having two VL 
derivatives on (f>, we only need 1 + e o/ them. We then go to the proof of Proposition 16 and note that it can 
be reproved assuming only that (p is in H^~^^ initially with suitable decay. 

Now Theorem 10 follows directly from Propositions 16 and 14. 

6 Estimates near the Event Horizon 

In this section, we will use the vector field Y to prove that any decay estimates that can be proved on a 
suitable compact set holds also along the horizon. We will also show that these estimates control enough 
derivatives to give pointwise decay estimate. 

Proposition 17. Suppose J^2((i.2)ro)*-r* {(f,) n'^ dV ok < Bt'" for all t, for some a>0. Then 

I Jl W <u +r')dVol^v=i(t,+r5)} <c(b+ [ {<t>) nldVolt^ 

Remark 9. The reader should think of B as some energy quantity of the initial data. For example, as we 
will show later, the hypothesis of this Proposition holds for B — C Ei {d^^Of'cf)). 

m+k<l 
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Proof. Apply Proposition 1 for Y, on the region 7?- = {5 (fi + rg) < w < | (f* + Tq) , t > fi} as in the figure, 
we get 



+ / Jl {4>) n>'^dVol{u=oo} + [ {4>) dVol 




Figure 3: The region TZ 



We spht up the integrals into r < ro and r > ro parts. 
Notice that the domain of integration of /r,>( 1 J}^ i4') n-^i 1 ,\dyoL ir^ , lies inside {r < rg}. More- 

over, we note that 1 ((^+^.j<„< 1 (■j^+^.-jj (0) n^dVo/{„=oo} > 0. Hence 



+ 

We estimate three terms on the right hand side. Notice that Y is constructed to be supported in {r < 
(1.2) ro}. 



/ Jl{4>)nldVolt, 

•> {\{ti+r'o)<v<\{t,+ro)} 
/■((l-2)ro)* 

<C Jl{<l>)n'i^dVok, 

J rX 



'0 

For the second and third term, we first use the compact support of Y and then apply the conservation law 
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associated to the Killing vector T. 

tY 



■/{r>ro}n{«>i(ti-r*)} 2(,*-+'^oJ ^ 2^ toJS 

<C [ J^{4>)n^^dVok, 

y{r5<r*<2((1.2)ro)*-r*} 



<CBtZ' 
<CBt^' 



[ \K^ {(j)) \dVol 

JTlr\{r>ro} 

.2)ro)* 

\K^ {(j)) \dVol 





r((l-2)ro)* 








''•0 




H(1.2)ro)* 


Itl ^ 






'"0 






<CB [ 


t-"dt 


Jti 





Jj {(j)) n'^dVoltdt 



<CB {U-ti)t^", 

since a > 0. 

Write / (t) = (0) <^^^^.^dT/oZ^„^i(,^^.^j. Then we have 

/ /(T)rfr<C7(/(ii) + Bmax{i*-ti,l}tr"). 

Jti 

We take C to be fixed from this point on. We clearly can assume the C > 1. 

^From this, we will prove the Proposition by a bootstrap argument. Assume / (t) < At~" for some large A 
that is to be determined. Wc want to show that / (t) < ■§t~°'. 

Let ti = — 8C^. Since wc arc only concerned with i« large, we assume without loss of generality that 
U>&[l-2-^^ so that t» < 2iii. Then 

r fiT)dT<C{Aq" + 8C^Bt^") 
<2C {A + 8C^B) t-"^. 

There exists i with — 8C^ <t<t>, such that 

t. 



{A + 8C^ 
- AC 



Now we let ti = t. Notice that U < 2at. Then 



/ {U) + j^' f (r) dT <C (/ (i) + 8C^Br^) 



<^C" + 2C2St7« + IQC^Btl 
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if A > 72C^B. 

Of course to have / {t) < At~°' for all t, we also need it to hold initially, i.e., A> f (to)- Therefore, we have 
/ i^) <U +,.)rf^o^{.= i(*.+.o)} <c(b+ f j;^ (0) nldVok^ t-", 

./{r<ro} 2V*0/ V / \ J J 

where C is a universal constant different from the one above. □ 

Using Proposition 17, we claim that a similar estimate holds on t-slices. 
Proposition 18. Suwose J^2((i-2)ro)*-ro* jt (^^j^^n'^dVolt < Bt'" for all t, for some a>l. Then 

[ Jj {ct>) nt^dVolr <c(b+ ( {ct>) n'^JVok,) y-", 

J{u,<-y<D,+l} \ J / 

for v*> 1. 

Proof. We prove this using the conservation law for Y on the region TZ = {i>, < v < + 2t!« — < t < r}. 
/ {<t>) K,+idVol{y=y,+i} 

J{t>2v,-r^} 

+ [ {(l))n'^dVolr+ [ K^{(t>)dVol 

J {v,<v<v,+i} Jn 

J{t>2v^-r^} J{v,<v<v^ + 1} 

We split up the integrals into r < ro and r > ro parts. 

Notice that the domain of integration of /|j>2?; -r'} ''T'l;^,dVoly^ lies inside {r < ro}. Notice also that 

/ Jj (<^)<.+idFoZ{„=„,+i} + / K''{cP)dVol>0. 

./{f>2'u»-r^}n{r<ro} Jnn{r<ro} 



Hence 



J^{cj>)n';dVolr 

{ti*<ti<?;* + l} 



< / < (<^) n^dVol,, + [ Jj (,^) n^2v,-r'dV0l^t^2v,-rS} 

J{t>2v,-r*} J{v,<v<v, + 1} 

+ [ Jl (0) <+WV^o/{„^„.+i} + / |if ^ (</.) \dVol. 

J{t>2v,-r5}n{r>ro} J7?.n{r>ro} 

We show that each term has the correct bound. The first term is bounded using Proposition 17, 

/ Jl {(fX^dVol,, 

J{t>2v,-r^} 



{r<ro} 



Jl{<i>)n'i^dVol,, 



<c[b + I (</.) nldVok,^ {2v. - rS)- 
<c{b + j (<^) nldVolt^j 
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The second term is controlled by assumption 



J{v,<v<v, + 1} 
J {v,<v<v, + l} 



<CB {2v,-rl)-" 
<CBv-°'. 

The last two terms are bounded by noting that Y is supported in r < (1.2)ro. The details are identical to 
the proof of Proposition 17. Therefore, 

/ (0) n>^dVolrC {b+ [ {ct>) n^t^dVok)^ t;-". 

J{v,<v<v, + l} \ J / 



□ 



This, and Sobolcv embedding, is sufficient to show pointwisc decay of the derivatives of along the 
horizon. We show further that if on a compact set, we have both energy decay and decay, then we have 
pointwise decay along the event horizon. More precisely, we have 

Proposition 19. There exist f very close to 2M such that if 

1 3-m .2((1.2)r-o)*-r* 

E E / i^l i^^^'^) < + '^') dVok < Bt- 

for all t, for some a > 0, then 

(1 3-m \ 
5 + E E / ■^/^ {9r^'<P) <dVok, v:'^, 
m=0 fe=0 / 

for v-t > 1, r < f . 

Proof. We first take f to be small enough to apply Y, i.e., f < tq. The exact condition on f will be determined 
later. 

For decay of (j) {v^ , r) , we want to show that on any time-slice, say t = t, 

For decay of dr* ^ (w* , r) , we want to show that on any time-slice, say t = t, 

/ {\f''dr^cj>f + \f''d^,<pf) dAdr; < Cv-'^. 

k^oJ{v,<v<v, + l} 

Proposition 18 gives 

/ ({dr,<l>f + \f(f>f) dAdr; <c(b + [j^ {<!>) n'^^dVok,) y-". 

J{v,<v<v. + 1} ^ \ J J 

After commuting with an appropriate number of Proposition 18 gives 

/ i\1Idr,4>? + mdr'4>\^) dAdr% <C[B + Y, [Jl^ (f^V) ) 

J{v,<v<v, + l} \ ^.^j^J J 
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After commuting with dt and using the equation, Proposition 18 gives 

/ ({dUf + \fdU\' + \rdU\') dAdr; 

J{v*<v<v^ + 1} ^ ^ 
/ 1 3— m „ \ 

\ ™ — n ; — n / 



m=0 k=0 



Therefore, it remains to show 



/ (j)'^dAdr* < Cv-". 

J {v,<v<v, + l} 



We rewrite 

p p2v,-t+2 

/ <l)'^dAdr; = / (j)'^dAdr*. 

J {v,<v<v,+l} J2t)*-T 

To achieve decay, we integrate in the w-direction and use the estimates we have on the compact set. 

2t),-T+2 



dAdr* 

r*+2 i-T i-2vt-t+2 



r-ZV 

J2v, 

< \ <^2^Adr|j=2«^-f}+ / / (pidu^) dAdr* dt 

Jf* * J2v,-f* J2v,-t 

ff*+2 PT p2v,-t+2 ,^ ,\2 „T p2v,-t+2 

< (t>'^dAdr\^^^^^_f ,^+ / dAdr*dt+ / / cf^ {I - n) dAdr* dt 



l2v,-f* J2vt-t ^ M J2v,-f* J2vt-t 

c2v,-t+2 pT p2v,-t+2 



pT p2v,-t+2 PT p2v,-t+2 

<B {2v^ - + / / K^{(j))dVol+ / / <I)'^{1- iJ,)dAdr*dt. 

Using the conservation law for Y, and controlhng all the terms on the region {r > f} with the assumption, 
we have 



PT p2v,-t+2 

(0) dVol 

J2v,-f* J2v,-t 
{r<r} 



where in the last step we have used Proposition 17. 
Therefore, 

p2v,-t+2 f [ \ p2v,-t+2 

/ (p'^dAdr; <C[B+ / ((^) n'^JVolto ] v'" + / (1 - dAdr*dt. 

The decay from the last term comes from the exponentially decaying (towards r* = — oo) factor (1 — /x). To 
use this decay, we use a bootstrap argument. Assume the decay Jy'I^~^^ (/f'^dAdr'^^.^^y < Av~", independent 

of t (Note that we can do this initially (in v) independent of t because after we fix f , the region of integration 
is a bounded set of the manifold. The apparent infiniteness is just an artifact of the choice of coordinates). 
We want to show that J^'_^~^^ (j)'^dAdr*^^^^y < ^v^'\ 
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^From the above, and using that (1 — /x) < Ce'^'' we have 

i-2vt-T+2 



'dAdrl 



J2v,- 

- 2 

if we choose A > 4C + / nt^dVok^) and P < -2 + i log □ 

7 Proof of the Main Theorem 

7.1 Improved Decay for 

To prove Main Theorem 1, we proceed in two steps. First, we show that for every there exist ti < t*, 
t such that a weighted L^-norm of (j) on the shce {t = ti} has the desired decay of t~^~^^- We then use 
the estimates for ip to upgrade this to decay estimates for a weighted L^-norm of (f) on the slice {t = t*}. 
We first set up some notation. Fix rj,r5. These are the rl and r'2 in the statement of Main Theorem 1. 
In other words, we would like to prove a decay estimate on the fixed compact region rf < r < Let 
> 2 (|rj| + jr^l) be the time slice on which we want to show the decay estimate. Let i = {1A)~^ and 
V ^{i<t<t^,rl-t^+t<r* <r*2+t^- t}. 

Proposition 20. There exist a ti with t < ti < t^, such that 

[ {(f + {dr, <pf) (1 - dAdr* < Ct-^Eo (0) . 

Jvn{t=ti} ^ ' 

t, 



u — const. Z. ^ X V — const. 



t^i^ (1.1)-% 



Proof. By Theorem 5, 



Figure 4: The region V 

j + {dr- (l>f) (1 - m)^ dAdr*dt 

<C I \K^'^^ (</-) I + V K"""^"' (00 dVol 

•It" I 
<Ci-^Eo (<^) 

<Ct:^Eo (</)) . 



Now take ti such that 



/ + {dr-(l>f) (1 - iJ,f dAdr* = inf / U'^ + {dr*(f>f) r'^ (1 - /if dAdr*, 

Jrn{t=ti} ^ ' i<t<t, J'pn{t=t} ^ ' 
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which exists since we are taking the infimum over a compact interval, and note that 

inf / ((j)^ + {dr'(l)f) r-^ (1 - dAdr* < (U - t)^^ [ {(f^ + (a^*</>)^') r^^ (i _ dAdr*dt 

i<t<t, J-pn{t=t} ^ ^ J-P V y 

<Ct-i J (j)'^ + {dr*(l)f') (1 - nf dAdr*dt 
<Ct-^Eo (</.) . 



□ 



To upgrade this to an estimate for a generic t, we make two observations about S. Firstly, S is timelike 
away from the event horizon. Secondly, S has a weight ~ t. We can therefore integrate from the "good slice" 
t = ti to the slice t = t^, and get the same decay estimate. This is done using integration by parts in the 
following Proposition. We prove a more general form but the reader should keep in mind that we will use 
f = ct>^ + {dr^cl>f,g = r-^ 

Proposition 21. Let f = f {r* ,t,ui & S^), g = g{r*), V = {ti < t < U,rl - U+t <r* <r^+U-t}. 
Then 

U [ fgdAdr* + [ vfgdAdr* + [ ufgdAdr* 

=h f fgdAdr* + 2 [ fgdAdtdr* + [ r* fg'dAdtdr* + [ {S f) gdAdtdr* . 

Jvn{t=ti} Jv Jv Jv 

Proof. We change to the variables u, v and integrate by parts, 

/ v{dyf) gdAdtdr* 
Jv 

= / v{dyf) gdAdvdu + / v (dyf) gdAdvdu 

J ^(t.-r^) Jti-u J ^(2ti-t,-r^) Jti-u 

= - I fgdAdtdr* - [ v fdr* gdAdtdr* + [ vfgdAdr* + [ vfgdAdr* 

Jv Jv Jvn{t=t,} Jvn{v=t,+r^} 

- / vfgdAdr*, 
Jvn{t=ti} 



[ u{duf) gdAdtdr* 
Jv 

= / / u{duf)gdAdudv+ / / 

J ^[iti-tt+rl) J ti-v J\[t,+rl) Jti-v 



u {duf) gdAdudv 



= - fgdAdtdr* + / uf dr' gdAdtdr* + / ufgdAdr* + / ufgdAdr* 

Jv Jv Jvn{t=t,} Jvn{u=t,-ri} 

- / ufgdAdr*. 

Jvn{t=ti} 

The proposition is proved by adding these two equations. □ 

To prove the main theorem, we use the above identity using f — (j)'^ + {dr*(j)f , g = (1 — jif. 
Wc notice that since /, g > by definition and u,v > Q m V = {ti < t < t^,r\ — t^ + 1 < r* < r2 + 1^^ — t} . 
Therefore, 

/ vfgdAdr* + / ufgdAdr* > 0. 

Jvn{v=i{t,+ri)} Jvn{u=^{t,-r-)} 
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Thus Proposition 21 would imply 

U f ((f^ + {dr*(l)f) (1 - pif dAdr* 

Jvn{t=t.} ^ ' 

^ (dr.^f\r-'^ {\- [if dAdr* ^2 [ ((l>^ + {dr*(f>f)r-^{l-fifdAdr*dt 
\r* (</.2 + {dr^(t)f^ (1 _ ^f^' \dtdr* + j + \dr.ipdr-'(t)\) r-"^ (1 - \if dAdtdr* 

<Ct-^Eo{(f>)+ I (j/j^ + {dr*^/jf^r-^ {1- fif dAdtdr* + f (j)^ + {dr*(f>f^ r-^ {1 - fif dAdtdr* 



<Ct:^+'Ei (0) , 



where we have used Proposition 20 at the second to last step and Theorems 5 and 10 at the last step. 
Therefore, 



' ((l>{t*f + {dr-(j){t.,)f'^ dAdr* < Ct-^+^Ei (</>) 



Since dt, are Killing, it follows immediately that 
1 



J2 r {{drdi.cj)iu)f + \f%.^it,) p) dAdr* < ct-^+'E, (arf^'» , 



1 = 



for any k, m. 

Using the equation Ogip — 0, we get 



f ' {W^cb{t,)f dVol < Ct-^+^ ^1 i^^^ 



The pointwise decay part of Main Theorem 1 follows form the standard Sobolev Embedding Theorem and 
Proposition 19. 

7.2 Improved Decay for dtcf) 

To estimate the time derivatives of 0, we follow an idea of Klainerman-Sideris [15j . The key observation is 
that the first derivatives oidt4> ^-re controlled with a weight of j^—^ by a linear combination of first derivatives 
of (j) and ip. This extra weight would give extra decay to dt(j). 

Proposition 22. Suppose t + r* > max{|, ^-^}- (This is true for example when r* is bounded below and t 
is sufficiently large.) 

1. \ {t- r*) 9201 < c {\dti>\ + \dr^i>\ + \dM + \dr-cp\ + (1 - /i) |r* 11^01) , 

2. |(i - r*)dr.dtC^\ <C{\dti'\ + \dr-.i^\ + \dt(j)\ + \dr>^\ + (1 - //) 1^11^01) , 

3. \t (1 - /i) < c((i - m) ly^l + \dr'n<j,\. 

Proof. Define A^^ = (1 - ^i)^^ dr-(l)+ f 9^*0 + ^0. Then Dg(j> ^ reads (1 - ^) d'^(j> = Ag(j). 
Recall that 

ip — tdtcf) + r*dr'(j). 

Therefore, 

dt^P-dt^^td^^ + r*dr^dtq}, 
dr* ■0 — dr* (j) = r* d'^,(f) + tdr* dtcj). 
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Hence, 

t [dtip - dt(f)) - r* {dr*i> - dr'cj)) 

= {t' - {r^f) dU + {r^f ((1 - m) A,<^ - e<^) 

= [t' - (r*f) dU + {r*f + (1 - ^) . 

Therefore, by re- arranging and dividing by {t + r*), 
\{t-r*)dU\ 

{t {dti^ - dt<f>) - r* {dr,i> - dr,<f>) - {r*f f + (1 - H) 



't + r 

<C{\dti,\ + + \dt4>\ + \dr'<l>\ + (1 - |r*||^<A|) . 

We have thus proved 1. On the other hand, using again the above equaUty, we also have 
{t-r*)dr.dtcl) 
= -dtip + dr*ip + dt(j> -dr*(j) + tdffj) - r*d%(l) 

= - + dr'i} + dt(l> -dr-4>+{t- r*) d^(f) + r* ((1 - Ag(f) - d^,(t)) 

2r* (1 — ij) 

= -dt^ + dr'i, + dt^ - dr.^+ {t - r*) d^(f) + ^ !^dr'(l)+ (1 - m) \r*\H- 

r 

This, together with 1, imphes 2. 

The proof of 3 is more direct. Using the definition of S, and that fl is independent of t and r* 

fl^ = r*dr'fl<f) + tdtfl(f>. 
Thus, by noting that O and differ only by constant, 

\t{i-n)fdtci^\ 
<\{^-^^)fi;\ + \{l-^l)-^r,^^(l>\ 

r 

<C{{l-fl)\f^\ + \drM<j>\). 

Corollary 23. 

Jj (M) nldVok, < Ct-^+' ( E E ^0 (aroV) + (<^) ] , 
for all c <1 and r. In particular, f can be chosen as that given by Proposition 19. 

Proof. We can consider large enough so that firstly, the assumption of of Proposition 22 holds and secondly, 
on the domain of integration, (f* — r*) t*. 

j"^' Jl{dt(t^)nldVolt^ 
['^' ({dr^dt<l>f + [dUf + (1 - m) \fdt<l)f) dVolt, 
<Ct-^ C' {{dti'f + {dr'^f + (1 - + {dt<t>f + {dr'<l>f + W^dt<t>\' + dVok, 

<CtZ^ / J"^ (V') + (<?i) + J"^ + J'^ {dtflcP) dVok, . 

J r* 

The corollary follows from Theorem 5 and 10. □ 



□ 
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Corollary 24. 

(2 4-m 2 2-m \ 

m=0 fe=0 TO=0 k=0 J 

Proof. We prove a Sobolev-type inequality. We first work on R^. We claim that for u G C° 

IHk~(«3)<qHl|,(^3)IHi(K3)- 

We give a simple proof using Littlewood-Paley theory. Let N & 2^ he a. dyadic number, x (0 be a radial 
cutoff function which is supported in {|^| < 2} and is identically 1 in {|^| < 1}. Define the Littlewood-Paley 

operators Pn by Pnu = (x [n) - X ( w)) u. 

Since the inequality claimed is invariant under scaling u {x) Xu {x) and u{x) ^ u (Xx) , we can assume 
that ||w||iji(R3) = ll''^llij2(R3) = 1. Then, by Bernstein inequality, 

\\PNu\\mm) < mm{CN-\CN-^}. 
Therefore, by Bernstein inequality again, 

iiwiiL=o(R3)<cE^^ii^iv«iu=(R3)<c[ E ^"'+ E 

N \N>No N<No I 



We note that a variant of this is true. We have for u G C', 



This is true because one can extend u into Bj. (0) without increasing the or norm. 

We now apply this to a cutoff version of 9^0. 

Let X — "I Q ja-j > 1 1 ' t = t^, \ct (f> = X for r > f, where f is as in Proposition 19, 

The norm is controlled with Corollary 23 and Theorem 5.4. 

\\dt4>\\HHm\Bm) 1/ (^*'^) <dVolt, + Ct:^ [ ' {dt<f>f (1 - dAdr* 

\m=0 fe=0 / 



The H norm can be controlled similarly once we note that using the equation, we have for r > r, 

\d^r'dtcl^\ <C {\dU\ + \\dr'dtcl>\ + \Tdtcl)\^ 
<C{\dU\ + \dr'dtcl)\ + \f^dt<i>\). 

Therefore, for t = t^, 

ll^«?^llij2(M3\B^(0)) 

0.55t» 

+ Ct:^ £ ' [{dfc^' + (dtQcpf + {dt<pf) (1 - m) dAdr* 

<ct-'+' f E E^o {dri^'ct>) + E E^^i {dri^'<i>)] ■ 



\m=0 k=0 m=0 k=0 
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Therefore, 

ll^t?^llL~({f*<r*<^}) <Pt0||L~(R3\B^(O)) 

~ 1 ~ i 

<<^Pt<^ll^l(R3\B,(0))Pt<^ll^2(R3\B,(0)) 

/ 2 2 11 ^ 



^m=0 fe=0 m=0 fe=0 



In particular, for t sufficiently large, this L°° estimate holds on sets of compact r*. Noting that the L° 
norm controls the L'^ norm on compact sets, we have 



L 



((1.2)ro)* / 2 2 11 

{dt<pfdvoit, < ct:^+' E E^o (^r^^V) + E E^i (sr^^V) 

\m=0 k=0 m=0 fc=0 / 



We also have, by Corollary 23, 

1 3-m ^((l.2)ro) 



E E / j^idrn'iM)) n'idVok, 

n 7 n J f 



m=0 fe=0 

(2 4-Tn 1 4-m \ 

^5:^0 (srov) + E E ^1 (^r^^'-^) ■ 
m=0 fe=0 m=0 fe=0 / 

The corollary then follows from the Sobolev Embedding Theorem and Proposition 19. □ 
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